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spacings, and obtain the preliminary result (/a: = 1.1872(41) in the continuum limit. 
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1. Introduction 

A great deal of interesting physics is accessible from the detailed study of the decay constants 
and masses of the light pseudoscalar mesons. For example, a determination of the decay constant 
ratio /k/Iti provides a competitive determination of the CKM matrix element \Vus\, but the lattice 
errors limit the precision of the result [1]. Similar fits can give ratios of light-quark masses, and 
with nonperturbative or perturbative renormalization, absolute values of the masses themselves. 
In addition, chiral fits make possible the determination of important low-energy constants. Here 
we present a preliminary study of light pseudoscalar meson masses and decay constants using the 
highly improved staggered quark (HISQ) action [2]. The MILC collaboration has been generating 
HISQ configurations [3, 4], which are expected to improve upon results from the "asqtad" staggered 
quark action. We attempt to fit the light pseudoscalar data with SU(3) staggered chiral perturbation 
theory (ChPT). Since data exists for light- and strange-quark masses that ai^e very nearly physical, 
we also perform simple linear interpolations of the data to the actual physical point and obtain a 
preliminary result for {fK/ fnY^^^^ ■ 

2. Ensembles and Correlators 

Several ensembles have been generated by the MILC collaboration using the HISQ action, 
with 2-1-1-1-1 dynamical flavors and lattice spacings ranging from approximately 0.06 fm - 0.15 
fm. The ensembles used in this report are found in Table 1. We use a notation in which mi is the 
common mass of the up and down sea quark, the mass of the strange sea quark, and the mass 
of the charm sea quark. Note that two ensembles have nearly physical m/, m.,, and m^.. Valence 
quarks are called x and y, and their masses are denoted by and lUy respectively. Pseudoscalar 
mesons with valence quarks x and y have mass ruxy and decay constant fxy. 

Quark propagators were generated with various sources and sinks, with 8-10 different 
valence-quark masses depending on the particular ensemble, and in the approximate range mf^^^ 
- ms'^-^\ where mf^-^^ and mf'''^' are the physical light- and strange-quark mass. (Propagators with 
valence-quark mass 0.9wic''^'' and mc'"" were also generated for heavy-light physics, but were not 
used in this report.) From these propagators, pseudoscalar meson correlators were calculated with 
random-wall and Coulomb-wall sources, and point sinks. The random-wall and Coulomb-wall 
correlators were fit simultaneously with independent amplitudes and a common mass. Only the 
amplitude from the random-wall correlators was used to obtain the decay constant. Four different 
time slices on each lattice were used for the source locations, with the locations of these slices 
shifting as we went through successive configurations. 

The fits started at a time tmin from the source con^esponding to ~ 2.4 fm in physical units 
for all of the lattice spacings. This choice was made by examining, on all of the ensembles, the 
dependence of fitted quantities on Fits including excited states and states of the opposite parity 
("alternating states") did not improve the quality of the fit, so they were not used. The exception was 
one of the ensembles with a nearly physical light-quai^k mass, which required adding an alternating 
state. 

We attempt to account for autocorrelations in the covariance matrix of the fitted masses and 
decay constants by using the following scaled covariance matrix in our chiral fits: 

r^f^Qj^, where rf^ = (2.1) 
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a (fm) 


ami 


arris 


arric 


Nlats 


0.15 


0.013 


0.065 


0.838 


1020 


0.15 


0.0064 


0.064 


0.828 


1000 


0.12 


0.0102 


0.0509 


0.635 


1040 


0.12 


0.00507 


0.0507 


0.628 


1020, 1000, 340* 


0.12 


0.00507 


0.0304 


0.628 


1020 


0.12 


0.00184 


0.0507 


0.628 


500** 


0.09 


0.0074 


0.037 


0.44 


1011 


0.09 


0.00363 


0.0363 


0.43 


633 


0.09 


0.0012 


0.0363 


0.432 


248** 


0.06 


0.0048 


0.024 


0.286 


503 


0.06 


0.0024 


0.024 


0.286 


280 



Table 1: HISQ ensembles used in this report. A single star indicates three different ensembles that have all 
the same parameters but different spatial volumes (L/a = 24,32,40). Two stars indicate ensembles where 
all quark masses are nearly physical. Note that all strange and charm quark masses are tuned to be nearly 
physical, except for the third ensemble with a = 0.12 fm, where the strange-quark mass is taken to be 
approximately 3/5 the physical value. 

lb] 

Here C/y is an element of the original covariance matrix, and C-^ is the corresponding element 
of the covariance matrix after blocking the configurations with block size b. Note that the scaled 
covariance matrix only uses the diagonal elements of the blocked covariance matrix Q . This is 
necessary because we use blocks of size b = 10, and not all of the off-diagonal elements of c|y' can 
be determined with good accuracy. 

3. Scale Setting 

Before proceeding to do combined chiral fits of all the ensembles, it is necessary to determine 
accurately the relative scale between ensembles with different lattice spacings. For each lattice 
spacing we take the ensemble with mj = O.lnis, and nig near the physical strange-quark mass. We 
look at masses nixx, and decay constants f^x, i.e., for pseudoscalars with two valence quarks of 
equal mass nix. Looking at the dimensionless ratio fxx/nixx as a function of rrix, we find the value 
of rux for which fxx/nixx = R, where R is some fiducial value. We take R = 0.3575, which is the 
value of fxx/^xx at nix = 0.4mf''^\ nii = O.lmf''", and nis = mf^-^'\ as determined from a previous 
analysis using the asqtad action [5]. We give to this value of Jiix the name "nip^s", or just nip for 
short, and the values of fxx and nixx at this point are called "/pp" and "nipp". Thus we expect to find 
that nip « 0.4m.v on each of the HISQ ensembles we are using to set the relative scale, and the extent 
to which this is not true depends on the accuracy of the tuning of rris to the physical strange-quark 
mass, an issue to which we will return later. (We ignore any possible mistuning of the charm mass 
as a small effect). 

In practice, we determine anip by looking at the three sets of valence masses anix = O.Sanis, 
OAarHs, and 0.6anis, and putting a parabola through m^xl fxx '^"^v- Once we have determined arup 
and afpp for a given lattice spacing, we can divide afxy and anixy by afpp for that lattice spacing, 
and similarly we can divide quark masses anix, aniy, ami, etc., by amp. This puts everything in "p4s 
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units". Quantities in p4s units are dimensionless quantities that are comparable between ensembles 
with different lattice spacings.^ With all quantities in p4s units, a combined chiral fit can be carried 
out. The advantage of our method is twofold. First, since the unphysical "p4s point" is well within 
the parameters of the data set, it can be reached by simple interpolation. Second, masses and decay 
constants at this point tend to have smaller errors than other quantities that could be used to set the 
scale. Our method is similar to the scahng trajectory method of [6]. 



4. SU(3) Staggered Chiral Fits 

We use SU(3), partially quenched, rooted, staggered ChPT to fit the pseudoscalar masses 
and decay constants as a function of valence- and sea-quark mass. The NLO expressions for the 
pseudoscalar mass and decay constant are given in Refs. [7] and [8], respectively. In our fits, 
we also include NNLO continuum terms from Ref . [9] , and add even higher-order analytic terms, 
which is similar to the previous asqtad analysis of Ref. [5]. In order to obtain a reasonable fit to 
ChPT, we find it necessary to impose a cutoff on the valence quark masses: nix+my <0.6{mi + ms). 
Also, our correlation matrix is close to singular, thus we smooth it by replacing all eigenvalues 
below 10^^ with their average. This is done as a quick fix in order to obtain prehminary results, 
and in future fits we will thin the data set instead. 

We show the results of the fit on a subset of the data in Figure 1 . The fit seems quite good to 
the eye, however when the covariances are included we see that the fit is not satisfactory: we get 
X^/d.o.f. = 1747/380 = 4.6. If we do separate fits for each lattice spacing we get ;^^/d.o.f. = 1.9, 
4.7, 2.6, and 0.7 for a = 0.15 fm, 0.12 fm, 0.09 fm, and 0.06 fm respectively. As discussed in 
Sec. 3, we assume that nis is equal to the physical strange-quark mass when we set the relative 
scale, which can be checked by seeing if nip/ms = 0.4. In fact, we find 0.37 < mp/ms < 0.43 for 
the different lattice spacings, indicating some mistuning of m^. However, this cannot completely 
explain the poor fit because the ;^^/d.o.f. is still large for fits to ensembles at single lattice spacings, 
where no scale setting is necessary. 

Another possible reason for the lack of a good fit is that quark masses are too large for SU(3) 
ChPT to work well, at least at the order to which we are working. In the asqtad analysis, this 
problem was dealt with by doing an initial fit with only lighter quark masses to fix LECs up to 
NNLO (and thereby controlling the extrapolation to mj''"'^), and then using the full set of quark 
masses to interpolate in the region of mf''^^ [5]. For the data set of this report (see Table 1), we only 
have one lighter-than-physical value of nis, so we cannot do something similar, at least in the sea 
sector. New ensembles are being generated that will make this possible in future work. Another 
possibility is to use SU(2) ChPT, in which the strange quark is considered heavy, to fit the data. 
Staggered ChPT expressions for masses and decay constants have now been worked out for the 
SU(2) case [10]. 



That is, if we knew fpp in MeV, then any meson mass or decay constant in p4s units from any lattice spacing would 
only need to be multiplied by this single number to put it in physical units. A similar comment applies for nip and the 
quark masses. 
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• a = 0.15fm 

T a = 0.12fm 

« a = 0.088 fm 

■ a = 0.058 fm 



m /m 

X p 



(a) Pseudoscalar decay constant as a function of = 



B 4.75 



4.25 



• 


a 


= 0.15 fm 


T 


a 


= 0.12fm 




a 


= 0.088 fm 




a 


= 0.058 fm 




(b) Pseudoscalar mass squared over the sum of the va- 
lence quark masses as a function of ra, — niy. 



Figure 1: Plots of the pseudoscalar mass and decay constant, and the fit to the data using SU(3) ChPT. We 
show only pions in these plots, i.e. m^ = my which is shown on the x-axis. Data on the plots are for all of the 
different lattice spacings, (shown in different colors), but only for ensembles with m; — Q.lm^. 



5. Linear Interpolations Near the Physical Point 

In Table 1 we see that there are two ensembles with a nearly physical light sea-quark mass. 
Since we also have valence quarks with nearly physical light- and strange-quai^k masses, we can 
do a linear interpolation or extrapolation as necessary to the exact physical point. Figure 2 shows 
plots of mly and fxy vs. wj^ + niy for one of the ensembles with nearly physical sea-quark masses. 
The plot of is remarkably linear. The plot of /vy is linear for small quark masses, but when one 
of the valence-quark masses niy is heavier, the slope starts to depend on niy. The situation is similar 
for the other ensemble, not shown. 

The observations of the previous paragraph suggest the following analytic fit form for a single 
ensemble: 

2 j {amj:f{mjc,my) =Ai+Bia{mjc + my-2mi), {mx,my) neai {mi,mi), 
I [amK) {mj^-,my) =A2+B2a{mx + my-mi-ms), {mx,my) near {mi,ms), 

afn [ntx , my )=C\+Dia{mx + my - Imi ) , {m^, my ) near (m/ , m/ ) , ^^^^ 

afKimj(,my) = C2+D2a{mx - mi) +E2a{my - m,.), {m^.my) near {mi,ms). 

We determine the values of the unknown coefficients A, - using data for am^y and afxy at 
{mx,my) = {mi, mi), {mi,ms), and the nearest neighboring points (m^^my) = {mi,m^^^), {m^^\ms), 
and {mi,m^^^), where m'^' is near m; and ni^^) is near m.,. For the ensemble with a = 0.09 fm, we 
have ami = 0.0012, am, = 0.0363, am^i' = 0.0024, and ^ o.029. 

Keeping the sea-quark masses fixed by continuing to work within each ensemble, we then 
wish to determine amf'^^^ and am^'^^\ where mj'''-^* and are the quark masses that correspond 
to the physical point. Here we ignore the sea-quark masses entirely which introduces only a small 
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Physical Pion Physical Kaon 



-I 



I 



I 



I 
I 



0.04 
a(m +m ) 







^0.0012 


• 


^ 0.0024 


9 am 


^ 0.00363 


a am 


- 0.00545 


• ^y 


^ 0.00726 


a am 

J 


= 0.0109 


■ am 


= 0.0145 


■ am 

y 


= 0.0218 


• ^y 


= 0.029 


a am 


^ 0.0363 



0.08 



0.06 p 
0.05 - 
0.04 - 
gS-O.OJ- 
0.02 - 
0.01 - 



Physical Pion 



Physical Kaon 



▲ Unitary Pion and Kaon 



0.02 0.03 
a(m +m ) 



(a) Plot of /„ vs. nix + ™v Points corresponding to 
different fixed values of amy are indicated with different 
colors. 



(b) Plot of m^y vs. nix + iny. The unitary points for this 
ensemble are indicated in red. 



Figure 2: Plots of w^,, and f„ vs. mx + my for the ensemble with a = 0.09 fm, ami = 0.0012, ams = 0.0363, 
and amc = 0.432. The values of m^ + my corresponding to the physical pion and kaon mass are shown as 
blue and green dashed lines respectively. These were determined using Eq. 5.3. 



paitial quenching eiTor because they are akeady close to physical. We set 

r I phys phvs\ nexp / phys phvs\ exp 



(5.3) 



where , nij^ , and are the experimental values. We can solve the left-hand equation for 
the unknown value amf^^\ and then solve the right-hand equation to determine the unknown value 
ams'^^\ This allows us to calculate the dimensionless ratio 

J afji ( m j''™' , mf^-^^ ) 

Results for {fK/fnY'^^'^ are shown in Table 2 for the two different ensembles with nearly 
physical sea-quark masses. We also show a continuum extrapolation of (/^//jt)^''"* by assuming 
that it is linear in for non-zero lattice spacing. To estimate the effect of partial quenching due 
to the fact that we kept sea-quark masses fixed during the analysis, we perform a similar analysis 
with unitary points only, from sever al of the ensembles with a = 0.12 fm.^ We find {fK/fnY''-' = 
1.1963(25), in agreement with the value of 1.1950(23) for the a = 0.12 fm ensemble from Table 
2. This gives us confidence that partial quenching is a small effect compared to the statistical 
uncertainty. 



6. Conclusion 

We performed fits using SU(3) staggered ChPT to pseudoscalar masses and decay constants 
on HISQ lattices with several different lattice spacings. We were unable to obtain a satisfactory 

^The reason this is not done in the first place is that it is only possible for a = 0. 12 fm with the current ensembles, 
(there is no lighter-than-physical strange-quark mass for the other lattice spacings), and because the interpolations are 
between points that are significantly further apart. 
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a (fm) 


ami 


aAmi 


ams 






0.12 


0.00184 


0.00019 


0.0507 


0.0013 


1.1950(23) 


0.09 


0.00120 


0.00022 


0.0363 


-0.0003 


1.1913(16) 


Continuum extrap. 










1.1872(41) 



Table 2: Results of interpolations to the physical point for /a'/A- We perform a continuum extrapolation 
assuming that this quantity is linear in for non-zero lattice spacing. We also show the distance of the 
interpolation points from the base data point, i.e. Aw; = — m/ and Aotj = OTi'''"* — ?Mj. 

fit, which could be due in part to strange-quark mistuning, and also possibly to the sea and va- 
lence strange-quark masses being too large for ChPT to work well. New ensembles with lighter 
sea strange-quark masses will allow us to address these issues in future work by 1) allowing us to 
correct for strange-quark mistuning, and 2) allowing us to use only lighter sea- and valence-quark 
masses to fix LECs to NNLO, and then to add higher order terms to interpolate near the strange 
mass using the heavier masses, which is an approach similar to Ref. [5]. We will also try using 
SU(2) staggered ChPT, which treats the strange quark as heavy. We were able to obtain the prelim- 
mary result ifK/fnY''^' = 1.1872(41) using simple linear interpolations around data points with 
nearly physical quark masses. The en^or on our result is comparable to that on the average of lattice 
results: {fK/ fnY^^^^ = 1.193(5) [11, 12], however we have included only the statistical eiTor in our 
result. When the full fit to a much larger set of data points is performed, we expect to be able to 
obtain a rather precise result. 

Computations for this work were carried out with resources provided by the USQCD Collabo- 
ration, the ALCF, and NERSC, which are funded by the Office of Science of the U.S. Department of 
Energy; and with resources provided by the NCSA and NICS, SDSC, and TACC, which are funded 
through the National Science Foundation's Teragrid/XSEDE Program. This work was supported 
in part by the U.S Department of Energy and the U.S. National Science Foundation. 
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